MICZ-KEPLER = DYNAMICS ON THE CONE OVER THE 
ROTATION GROUP 



RICHARD MONTGOMERY 



Abstract. Wg show that the n-dimensional MICZ-Kepler system arises from 
symplectic reduction of a simple mechanical system on the cone over the ro- 
tation group SO{n). As a corollary we derive an elementary formula for its 
general solution. The punch-line of our computation is that the additional 
MICZ-Kepler |0p/r'^ type potential term is the rotational part of the cone's 
kinetic energy. 



1. Introduction 

The classical mechanical formulation of a Hydrogen atom is identical to the 
Kepler problem for a planet moving in the gravitational field of a massive Sun. The 
MICZ [Mclntosh-Cisneros-Zwanziger] -Kepler system is an integrable extension of 
the Kepler problem in which the charged proton at the origin of the Hydrogen atom 
is simultaneously a Dirac monopole. The original references are [Hill]- See |5J[5] for 
history, references, and the n-dimensional generalization. Our purpose is to show 
that the MICZ-Kepler system is the reduction of a natural mechanical system (no 
magnetic fields) on the cone over the rotation group. See theorem 1 below and eq 
([8]). We then use this realization to write down an explicit Lie-theoretic formula 
(eqini) for the system's general solution. 

Meng [3] formulated the MICZ-Kepler system as a system on what he called a 
'Sternberg phase space' -and what we will call the 'adjoint bundle phase space' - 
associated to a principal SO{n — 1) bundle over M" \ {0}. This phase space arises as 
the symplectic reduction of the cotangent bundles of the same principal bundle. See 
[3 [6l H] . It follows that there is a Hamiltonian system whose configuration space 
is Meng's principal bundle and whose reduction yields the MICZ-Kepler systems. 
Our object is to find this system, and then use it. The heart of the computation 
(noticing that the Hamiltonian eq([8|) becomes (|14p in adjoint bundle variables) is the 
observation that the angular part of the Kepler kinetic energy plus the quadratic 
color charge term of the MICZ-Keple Hamiltonian equal the kinetic energy for 
the bi-invariant metric on the full rotation group. (For the quantum-mechanical 
analogue of this computation compare the algebra around eq ((T4| to the algebra at 
the end of section 2 of [9]. ) 

2. MICZ-Kepler 

The MICZ -Kepler system of equations can be written as a mixed lst-2nd order 
system for a curve {q{t),(f){t)) in a vector bundle over R" \ {0}. Here q{t) denotes 
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the curve in R" and (j>{t) lies in the moving fiber over this curve. This fiber is the Lie 
algebra so{n — 1) of SO{n — 1). The bundle is called the adjoint bundle, and is an 
associated vector bundle to a certain principal bundle SO{n — 1) — > Q — > M" \ {0}, 
called the Dirac monopole and described in the next section. So, relative to a 
local trivialization of the bundle, takes values in the Lie algebra so{n — 1). 
The principal bundle is endowed with a canonical connection whose curvature is 
F, and which induces a covariant derivative D on the adjoint bundle. Then the 
MICZ-Kepler system is 

D(l)/dt = 0. 

Some clarifications are in order. The curvature is a two-form with values in the 
adjoint bundle. Consequently F(ci, ■) is a one-form with values in the adjoint bundle. 
The Killing form endows the adjoint bundle with a natural fiber inner product, so 
that (j) ■ F{,q,-) is a one- form on R". We turn this one- form into a vector using 
the standard flat metric. In coordinates then: 4> ■ F{q, •)■' = (f)aF^j{q)q'' . In a local 
trivialization the second equation reads D(j)/dt = d(j)/dt + [Aiq^, 0] where A is the 
connection one- form, an so(n — 1)- valued one- form , in this trivialization. 
The system has two basic conserved quantities: its energy 

(1) H^i,Mf^%-l. 

and the square norm of the adjoint variable, or so{n — 1)- 

Casimir= ||0(t)f . 

There are other conserved quantities, an angular momentum and a Runge-Lenz (or 
Laplace) vector, but we will not need them here. 

The adjoint bundle canonically fibers into adjoint-orbit fibers which are preserved 
by parallel transport. Consequently the fiber variable <j){t) stays on whichever 
adjoint orbit fiber it begins on at time 0. As a particular case, the zero-orbit is 
preserved and corresponds to setting (/) = in the equations, which reduces them to 
the standard equations of the Kepler problem. Meng [3] takes 4) to lie in an adjoint 
orbit of a particular type which he calls "magnetic" . We allow any 0, hence any 
adjoint orbit. 

3. The n-dimensional Dirac Monopole 

Wc describe the bundle- with-connection needed to define the MICZ-Kcplcr equa- 
tions of the previous section. We take our description from Meng, who called it the 
Dirac monopole, since that is what it is when n = 3. 

Consider the space Co of all orthogonal frames / = (/i, /2, ■ • ■ , /«) on M" nor- 
malized so that their lengths are all equal: = \fj \ all i, j. The map 

/ ^ /„, Co ^ M" \ {0}, 

gives Co the structure of a principal 5*0(71 — 1) bundle over R"\{0}. Write r = \ fn\- 
Then, we have an S'0(n)-equivariant diffeomorphism 

R+ X SO(n) -> Co;(r,,9) h-> {rgei,rge2, . ■ .rgcn) = (/i,.../„), 

where ei, . . . ,e„ is the standard basis of M". From this perspective, the bundle 
projection becomes {r,g) i~> rgcn- 
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Put spherical coordinates on M" so that R" \ {0} = R+ x 5"'"^ Recall that 
= SO{n)/SO{n - 1) where the SO{n - 1) C SO{n) is the stabilizer of 
e„. The restriction of the principal bundle Co ^ K" \ {0} to S'""^ C R" \ {0} 
defines the homogeneous principal bundle tt^^-i : SO{n) — >■ 5"^^. The Killing form 
on so{n) endows SO{n) with a bi-invariant metric cPsgo{n)i ^-nd relative to this 
metric, the orthogonal complement to the fibers of the projection 7r5n-i define an 
S'0(n)-equivariant connection for the homogeneous principal bundle. Extend this 
connection trivially in the radial (r) direction to arrive at the "Dirac connection" 
of Co — >■ K" \ {0}. We refer to the bundle Cg endowed with this connection as 'the 
monopole'. For coordinate expressions for the connection see [9l|4j[2l[3] . 

4. Adjoint and co-adjoint bundles. 

Suppose that G ^ Q — 5* is a principal C bundle. Then the adjoint bundle 
Ad{Q) := Q Xq g — > 5 is the vector bundle associated to Q via the adjoint action 
of C on g. This means we divide the product Q x g by the C-induced equivalence 
relation (q,^) ~ {qg,Adg-i^). Write equivalence classes [q,(,]. The fiber of the 
adjoint bundle is g. 

The co-adjoint bundle is defined similarly, using the co-adjoint action. If g 
is endowed with a bi-invariant inner product such as the Killing form on g = 
so(n — 1), then we get G-equivariant isomorphisms g = g* and hence we can identify 
the adjoint bundle with the co-adjoint bundle. This identification sends adjoint 
orbit fiber to the corresponding co-adjoint orbit fiber. We make this identification 
throughout the paper. Thus the variable in the MICZ-Kepler system is a 
section of the adjoint bundle Ad{Co) along the curve q{t) g K" \ {0} 

A connection on Q induces on any associated vector bundle, and so on the adjoint 
bundle. We can describe parallel transport in q{Q) along a curve c(i) G 5 as follows. 
Pick a point 0(0) = [(7(0),^] G q{Q)- Consider the horizontal lift q{t) G Q of c{t). 
Then = [q{t),^] G g(Q) is the parallel transport of (j){0) along c{t). The 2nd 
MICZ-Kepler equation states that is covariantly constant, and hence of the 
form just described. 

5. Metric cones and Riemannian submersions 

Let r — \fn\ — > in the construction of Gq = x SO{n) so that the whole 
of SO{n) is crunched to a point defined by r = 0. We we have formed the cone 
G — Cone{SO{n)) D Co over the full rotation group SO{n). The bundle projection 
Go R" \ {0} extends to the cone point r = 0, sending it to the origin. 

We now put a canonical metric on the cone C. 

Recall that the cone over a Riemannian manifold (X^d^sx) is given by the 
metric dr^ + r^d^sx on R+ x X where r G M+. As r the metric factor 
involving X shrinks to so it makes sense to identify {0} x X to a single point, 
which is the cone point. In this we get a metric on Cone{X) = [0, 00) x X/ {0} x X 
which is Riemannian away from the cone point. 

We form the metric cone G = Cone{SO{n)) by applying this cone construction 
to X = SO{n) endowed with a Killing induced bi-invariant metric d'^ssO(n). Such 
a metric is well-defined up to scale. That scale will be fixed by insisting that the 
bundle projection Gq M" \ {0} is a Riemannian submersion. 

Recall the notion of a Riemannian submersion t: :Y ^ S between Riemannian 
manifolds Y, S. Suppose that tt is a submersion. Take any point j/ G F and consider 
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the orthogonal complement at y to the fiber 7r^^(s) through y. Here s = 7r(y). 
Following the principal bundle language as above, call this orthogonal complement 
T-Ly C TyY the 'horizontal space' at y. The differential diTy of tt, restricted to 
the horizontal space, is necessarily a linear bijection onto the tangent space to 
s. Now the horizontal space inherits an inner product from Y. If this restricted 
differential is an isometry between inner product spaces for all y, then tt is said to 
be a Riemannian submersion. 

Consider the standard metric cPsg^-i on the unit sphere 5*"^^. The Killing 
scale for SO{n) is now fixed by insisting that the bundle projection SO{n) S""'^ 
be a Riemannian submersion. This scaling is the one for which the standard 
basis elements A ej of so{n) have length 1. Here a A ej is the skew-symmetric 
operator sending Ci to ej and ej to — e^, i ^ j. Let 9ij be the dual basis, viewed as 
left-invariant one-forms. Then (Psso{n) = ^{Oij)^so that 

(2) d^sc - dr2+r2S(%)2 

We verify that SO{n) — > 5*""^ and Co \ {0} are Riemannian submersions. 

The connection form for the Dirac monopole is the so{n — l)-valued one-form 

(3) A^T.i^j^nOije, Acj, ( on S'O(n) or Co). 

(Note we must use left invariant one forms, rather than right-invariant formsas 
can be seen by the fact that the connection is not C-invariant, but rather C- 
equivariant with G acting on the Lie lagebra by the adjoint action.) The horizontal 
distribution is defined by 9ij = 0,i,j ^ n. The vertical distribution is defined by 
Sin ~ 0,i — 1,2, . . . ,n — 1, together with dr = in the Cq case. Thus the Killing 
metric on SO{n) splits orthogonally relative to the horizontal- vertical splitting 

(4) d?SS0(n) — S(0i„)^ + Si<j<„(%)^ 

(5) = 7r*d^Ss„-l -I- d^S fiber 

The fact that exp(teiAen) has period 2tt shows that the metric scalings are correct 
for the Riemannian submersion: a 27r periodic horizontal geodesic in SO{n) maps 
to a great circle of circumference 27r in 5""^. 

Now, written out in spherical coordinates the metric on R" is c?Sg„ = dr^ + 
r^dsg„_i, which is to say that metrically speaking M" = Co7ie(S'"^^). The corre- 
sponding radial-horizontal spherical- vertical splitting of the metric on C is 

(6) d^sc = (rfr2+r2l](0„O')+7-'S,<j<„(%)2 

(7) = TT*ds'^r, + d'^S fiber 

which shows the projection C — > M" becomes a Riemannian submersion away from 
the cone point. 

6. Main result: MICZ-Kepler from a mechanical system on the cone. 

A natural mechanical system consists of a configuration space Q, endowed with 
a Riemannian metric dsg and a potential function V : Q ^ M.. This data defines 
a Hamiltonian on T*Q whose Hamiltonian H is kinetic plus potential: H = K + 
V, where the kinetic energy K is induced by the metric. In standard canonical 
coordinates {q,p) = {qi,p^) for T*Q we have K{q,p) = ^g^-' {q)PiPiPj if dsg = 
gijdq^dq-'. 
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Theorem 1. Take the canonical metric cone C = Cone{SO{n)) described in the 
previous section so as to get a Riemannian metric on Co = Cone{SO{n) \ {0}) — )■ 
R" \ {0}. Take potential function V ^ -\ : Cq where r : Cone{SO{n)) 

[0,oo) denote the cone's radial coordinate. Then the resulting natural mechanical 
system (Hamiltonian ^ below) on T*Co is invariant under the lifted action of 
SO{n), and so is also invariant by SO{n — 1). The symplectic reduction of this 
system by SO{n — 1) at any particular fi £ so{n — 1)* = so(n — 1) yields the 
generalized MICZ-Kepler system associated to the adjoint orbit through fi. 

We see from the expression ([2]) for the metric on the cone that the Hamiltonian 
of this theorem is 

(8) + 

where r,pr are canonical coordinates on T*R"*" and the arc the Lie-Poisson 
coordinates - linear coordinates on so{n)* - induced by the choice of basis A Cj 
for so(n). 

Corollary 1. Any solution to the generalized MICZ-Kepler system can 

be constructed as follows. Fix ^ G so{n),w € 5""^. Fix a solution r(t) to the 
1-dimensional Kepler problem: f ~ ^V^ii') where the effective potential is V^{r) — 

^7 + 72 with ~ l^p, and a solution u{t) to the ODE u ~ ^(fyj • Then 

(9) q{t) ^ r{t)exp{u{t)Cjw 

The adjoint bundle variable (f){t) is obtained by parallel translating an initial adjoint 
vector (t>{0) = [/, ^/(^)] along q{t), where f £ Co is any element projecting to q(0) 
and A is the connection one- form Q). 

Special Cases. 

1. Kepler. Take ^ horizontal over the initial (?(0), so that A(^) = and 
4>{t) = 0. For simplicity, take the initial ^(O) in the direction w ~ e„. Then 
horizontality implies ^ = Sw^ej A e„ = ij' A e„ is an infinitesimal rotation in the v, e„ 
plane. The solution q{t) then lies in this plane. Set 0{t) = u{t)\^\. Then {r{t),6{t)) 
is a solution to Kepler's equations expressed in polar coordinates. 

2. Magnetic cone. Meng takes his 0(0) to be of "magnetic type', which means, 
relative to a local trivialization, that (/)(0)^ = —fxld. In other words, up to scale 
(j){t) is an almost complex structure on K", compatible to the standard complex 
structure. Thus ^ satisfies = — /i^/d. Set J = ^//x so that J is an honest almost 
complex structure. Then exp{u^) ~ cos{pu)I -\- sin(nu)J . The solution q{t) lies on 
the two-plane spanned by q{0) and Jq{0). Indeed, it is another Keplerian conic on 
that plane, as Meng showed in [3]. 

3. Generic. Take ^ generic, meaning that it has [n/2] distinct nonzero eigen- 
values ±iujj linearly independent over the rationals. We can, by conjugating by a 
rotation, put ^ into the normal form TiUJje2j-i A e2j. Then 9 1— >■ exp{9S,) is a dense 
curve on a standard maximal torus in SO(n). For negative energy the correspond- 
ing one dimensional Kepler motion is periodic with period T , and without collision 
(since £, 0). We can arrange that 27r/T is rationally independent of the ujj. Then 
the corresponding solution curve q(t) forms a dense winding on a kind of "annular 
projection" r^in <r< r^ax to R" of a torus of dimension [n/2] + 1 . 
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7. Proof of theorem 1. 

Wc apply the general theory of reduction of cotangent bundles of a principal 
bundle. We first describe that general theory. See [S], particularly pp. 160-163, or 
the earlier references [Sj [U [7] for perhaps more leisurely descriptions. 

Let G ^ Q ^ S he a principal G-bundle. G acts on T*Q with G-equivariant 
momentum map J : T*Q g*. The quotient {T*Q)/G is naturally a Poisson 
manifold whose symplectic leaves are the symplectic reduced spaces J^^{0^)/G = 
J~^{li) /G ^ where = G ■ fi denotes the coadjoint orbit through /i G g*. The gen- 
eral theory proceeds by using a connection on Q to define a symplectic isomorphism 
with these reduced spaces. 

Differentiate the sequence of maps G ^ Q — S at fixed q E Q to obtain 

the sequence of linear maps g *- TQ TgS where s = i^yq). The sequence 

is exact: the image of aq equals the kernel of diTq. (This common image is called the 
vertical space at q.) Letting q vary parametrically we obtain the 'Atiyah sequence' 
(described in [1]) 

(10) Qxg^—^TQ—^TT*sTS 

which is an exact sequence of G-cquivariant vector bundles over Q. Dualizing the 
first map of (|10p . and composing with the projection yields the momentum map : 

J: T*Q^^QxQ* ^Q* . 

A connection ^ for Q 5* induces a G-invariant splitting of ([TU| and hence a 
G-equivariant isomorphism: 

TQ'=^tt%TS®{Qxq). 
Dualizing yields the G-equivariant isomorphism: 

(11) T*Q^Ti*sT*S®{QxQ*). 

Now G acts on the bundle Q x g by g{q,C) = {qg, Adg-i^) as per the equivalence 
relation used to define the adjoint bundle and thus the action of G on Q x g* is the 
one used to define the co-adjoint bundle. Forming the quotient by G we thus get 
the bundle isomorphism 

*A : {T*Q)/G ^ T*S ® Ad*{Q) 

over S which is our desired identification. We refer to the right hand side of this 
isomorphism as being "on the Adjoint bundle side" in what follows. 

From our factorization of J we see that under the isomorphism 'J'^ the reduced 
spaces J~^{n)/Gf_, = J~^{0^)/G become the submanifolds T*S® (0^)((3), which 
are the Adjoint bundle phase spaces and the Sternberg phase spaces of ([3]). If g is 
endowed with a bi-invariant Killing form as above, then the co-adjoint orbit bundle 
{Ofj){Q) is identified with a corresponding adjoint orbit. 

We need more detail regarding the isomorphism a and the Poisson brackets 
on the universal phase spaces in order to pull-back the cone Hamiltonian ([8]) and 
compute equations of motion. The connection defines horizontal lift operators hq : 
TgS — > TqQ, q Q, which are linear operators whose image is the horizontal space 
Kq = ker{A{q)) C TqQ of the splitting of TQ. The dual of hq is h* : T*Q T^S 
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and is one factor of the isomorphism ([TT]) . Write [q, P] for the equivalence class in 
T*Q/G of (g, P) e T*Q. Then 

^A{[q, P]) = (^(g), /i^P) e [g, J(g, P)] e T*S ® ^d*(Q). 

We describe "^a in coordinates. Let {x\g) e R'' x G be coordinates on Q in- 
duced by a local trivialization Qu = U x G oi Q S, together with coordinates 
on J7 C S". Then, over U wc have T*Q ^ T*U x T*G = T*U x G x g* with co- 
ordinates {x'\pi, g, £,a) where the coordinates ai'c Lic-Poisson linear coordinates 
on g*. relative to a basis Ca of g. Thus {T*Q)/G = T*U x g* with coordinates 
Xi,pi,^a- On the other hand, the same data {x\g) and basis Ca yield coordinates 
Xi, iTi:S,a for T*S ® ^c?* (Q). Relative to these two sets of coordinates the map ^ a 
is the minimal coupling procedure *A(a;,Pi,^) = {x^- ,pi- £,aA°;[x),£,a) = {x\7ri,£,a) 
where A{x) = 'SeaAf{x)dx is the connection one-form relative to the local trivial- 
ization and coordinates. The brackets on the Adjoint bundle side are {xi,T:j} = 
SjAUb} = -cUd and {tt.„tt,} = -^.J^g, {tt^^,} = AC = [A^,^]a■ Here i^^ 
is the expression for the curvature of A in this local trivialization. (Compare eqs 
(12.2) of [5] to (3.2) [3]. Note that in the triple of displayed equations immediately 
following (12.2) of [5] most terms should have a capital P immediately in front of 
them.) These agree with the brackets found in Meng for the case of so{n — 1). 

Now we recompute the Hamiltonian ([8]) on the Adjoint bundle side using '^a- 
The dual of the metric splitting ^ yields 

(12) Kc = ^[p^ + ^(S-LieL + ^S,<,<„4)] 

(13) = ljpl + ^ih*Ks^ + \cj^n 

where h* : T*SO{n) — !> T*S'"^^ is the connection induced dual of the horizontal lift. 
The fiber term Ei<j<„^,^j- corresponds, on the Adjoint bundle side, to the Casimir 
function of SO(n — 1), viewed as a function on the adjoint bundle. So the 
Hamiltonian, viewed on adjoint bundle side, reads 

H = \{pl + ^Ks^ + ^\<l^\^)-- 
It r r 

The sum of the first two terms ^(p^ -f -^Ksr^) is the usual kinetic energy ^E^^j^tt^ 
on K" written in spherical variables. Thus 

(14) i7= 

2 r 
which is the MICZ-Kcplcr Hamiltonian. 

8. Proof of the corollary. 

We compute the equations of motion on Go, using the expression (jH]) for the 
Hamiltonian and the fact that Vl = j|^||^ is a Casimir. Set /i^ = and V^{r) ~ 

-7 + 7^ Then Hamilton's equations on T*Go = T*M+ x T*SO{n) = M+ x M x 
SO{n) X so{n) are: 

r = {r, H} = pr 
p, = {p.,,H]^-V'^.{r) 
OH 1 , 
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The first pair of equations decouple from tlie second pair, and assert that (r^pr) 
evolves as per the one-dimensional radial Kepler equation with effective potential 
V^. The last equation asserts that ^ G so(n) is constant. The equation for g 
asserts that g{t) = goexp{u{t)^) where du/dt — Indeed, the solution to 

g = gS, through g = Id is the one-parameter subgroup exp{t^) and this flow is 
generated by the Hamiltonian ifi. We have scaled the Hamiltonian on SO{n) by the 
(time-dependent) factor and used left-invariance. To rewrite goexp{u{t)^) = 
exp{u{t)i)go we can set | = ^o^ffcT^- 

We now have the solution the Kepler equation on the cone: {r{t),g{t)) with 
g{t) — exp{u{t)(^)gQ. Recall the bundle projection is (r, g) i— > rgcn and use that any 
unit vector w can be written goCn to the expression for q(t) in the corollary. QED 

9. Other groups 

The tricks used here apply to any Lie group G in place of SO{n) provided that 
G is endowed with an faithful orthogonal representation on R" which is transitive 
on the unit sphere S*"^^. We get the theorem that the Kepler problem on Cone{G) 
is equivalent to the 'MICZ-Kepler-G ' problem whose 'color variables' lie in an 
adjoint orbit bundle for G over M" \ {0}. The standard families of such groups are 
the unitary groups U{n) on R^" ^ c" and 5p(n;H) on H" = R"". 
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